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Abstract 

We construct a set of exact ground states with a localized ferromagnetic domain 
wall and an extended spiral structure in a quasi-one-dimensional deformed flat- 
band Hubbard model. In the case of quarter filling, we show the uniqueness of the 
ground state with a fixed magnetization. The ground states with these structures 
are degenerate with the all-spin-up and all-spin-down states. This property of the 
degeneracy is the same as the domain wall solutions in the XXZ Heisenberg-Ising 
model. We derive a useful recursion relation for the normalization of the domain 
wall ground state. Using this recursion relation, we discuss the convergence of the 
ground state expectation values of arbitrary local operators in the infinite-volume 
limit. In the ground state of the infinite- volume system, the translational symmetry 
is spontaneously broken by this structure. We prove that the cluster property holds 
for the domain wall ground state and excited states. We also estimate bounds of the 
ground state expectation values of several observables, such as one- and two-point 
functions of spin and electron number density. 

keywords: ferromagnetic domain wall, spiral state, flat-band Hubbard model, ex- 
act solution, quantum effect, cluster property 
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1 Introduction 

Appearance of domain structure is known to be a universal phenomenon in ferromagnetic 
systems with an energy gap and with a finite correlation length. In realistic materials, the 
main driving force to produce the ferromagnetic domains is the magnetic dipole-dipole 
interaction, which is much smaller than the Coulomb interaction but long-ranged. In some 
examples we introduce here, however, stable domain structures appear without such long 
range interactions. Alcaraz, Salinas and Wreszinski construct a set of ground states with 
a domain structure in the ferromagnetic XXZ model with a critical boundary field in 
arbitrary dimensions for an arbitrary spin |T|. They have shown that the degeneracy of 
the ground states corresponding to the location of a domain wall center is identical to 
that of the ground states in the SU(2) invariant model. In their solution, the domain wall 
is localized at an arbitrary location with a finite width depending on the Ising anisotropy 
parameter {q+q~ 1 )/2 > 1 with a positive number q. After their construction of the domain 
wall ground state, lots of interesting properties of the solution have been discovered. 
Here, we focus only on one-dimensional models and remark several interesting results. 
Koma and Nachtergaele have proved the existence of the energy gap which is independent 
of the location of the domain wall J2] in the one-dimensional XXZ model with a spin 
S = 1/2. Matsui has obtained the complete set of pure ground states in this model in 
the infinite- volume limit [3]. Koma and Nachtergaele have also obtained the complete 
set of pure ground states with domain walls in the one- dimensional XXZ model with 
an arbitrary spin S in the infinite- volume limit J3]. Koma, Nachtergaele and Starr have 
shown a singular excitation localized near the domain wall with lower energy. Bach and 
Macris evaluate a spin one-point function in the domain wall state of the one-dimensional 
XXZ model by a rigorous perturbation method J3]. Datta and Kennedy discuss the 
stability of a domain wall in one-dimensional quantum spin models at zero temperature 
by another rigorous perturbation method Kirkwood-Thomas approach around the Ising 
model as an unperturbed model JH]. In the antiferromagnetic XXZ model, the exchange 
interaction destabilizes the domain structure, while in the ferromagnetic model it preserves 
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the domain structure at zero temperature. The role of the quantum effects should be 
studied more in many other models. 

In many-electron systems, stable ferromagnetism itself is less trivial than that in spin 
models. The ferromagnetism appears by purely quantum effects in many-electron sys- 
tems. There are no corresponding classical systems to the ferromagnetic many-electron 
systems. Some remarkable results for ferromagnetic ground states have recently been 
obtained in a class of Hubbard models. Mielke has shown that the ground states give 
saturated ferromagnetism in many-electron models on general line graphs with special 
properties [Zj. Independently, Tasaki has shown that ground states give the saturated 
ferromagnetism in another class of models jH] . These models are called the flat-band Hub- 
bard models, since the band structure for a single electron in the carefully constructed 
hopping Hamiltonian becomes flat. Mielke and Tasaki argue the sufficient conditions for 
the flat-band ferromagnetism in extended models [D]. Tasaki has proved that the sat- 
urated ferromagnetism is stable against a perturbation which bends the electron band 
Tanaka and Ueda have shown the stability of the saturated ferromagnetism in a 
more complicated two-dimensional model in Mielke's class [H]. The ferromagnetism in 
flat-band Hubbard models and their perturbed models is believed to be one universal 
nature of ferromagnetism in many-electron systems. In these models, the ferromagnetism 
appears in nonperturbative region in which the Coulomb repulsion between electrons is 
sufficiently strong. No perturbative expansion method can explain the ferromagnetism 
in these models. Only non-perturbative methods particularly devised can work for these 
problems. 

In this paper, we construct a set of exact ground states in a one-dimensional flat-band 
Hubbard model with coexisting domain wall and spiral structures. We deform a flat- 
band Hubbard model by introducing a complex anisotropy parameter q. The SU(2) spin 
rotation symmetry in the original flat-band model is reduced to U(l) symmetry in our 
deformed model. This anisotropy \q\ ^ 1 leads to a localized domain wall with a finite 
width and a complex q leads to an extended spiral state. We see that the domain wall 
ground state in our electron model has the same degeneracy as that in the XXZ model 
and the same localization property in a certain parameter regime. As in the XXZ model, 
the domain structure is characterized in terms of the local order parameter (S x ), which 
represents the third component of the localized spin at site x. This local order parameter 
takes the same sign within one domain. The domain wall center is a set of sites x defined 
by zeros of the local order parameter (S®) = 0. For complex q = \q\e>° with real B, 
our electron model differs from the XXZ quantum spin model which has no spiral ground 
state. The ground state in our model has an incommensurate spiral structure with a pitch 
angle 9. In the case of quarter filling, we show the uniqueness of the ground state with a 
fixed magnetization for an arbitrary q. 

To calculate the ground state expectation values of observables, we introduce a nor- 
malization function of the ground state. In particular in the one-dimensional model, we 
derive a recursion relation for the normalization function, which determines this function 
together with its boundary data. This recursion relation is a linear equation in the un- 
known function of a site with a function coefficient which depends on the distance from 
the domain wall center. This method is similar to the method for the normalization func- 
tion in the non-product domain wall ground states in the XXZ model obtained by Bolina, 
Contucci and Nachtergaele [T2*l ITS] . We treat this recursion relation by a kind of per- 
turbative expansion method in the function coefficient. This treatment of the recursion 
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Figure 1: The lattice A. White circles are sites in A D and black dots are sites in A'. 
Electrons at a site can hop to other site if the site is connected to the original site with a 
line. 



relation enables us to discuss the convergence of the ground state expectation values of 
local observables in the infinite- volume limit. Then, we prove that the cluster property 
holds for the domain wall ground state and excited states. This property implies that the 
obtained representation is healthy as a quantum system with infinite degrees of freedom. 
As discussed in |14j . we expect that quantum effects in the Hubbard model stabilize the 
domain wall as in the XXZ model. Finally, we show exponential decay of the truncated 
correlation functions, and evaluate their decay length. 

This paper is organized as follows. In section EJ we define a deformed flat-band model 
on a quasi-one-dimensional lattice. We give an explicit form of a Hamiltonian. In section 
EJ we construct ground states, and prove the uniqueness of the ground state with a fixed 
magnetization. We also introduce different basis of ground states whose state is super- 
posed over all magnetizations with certain weights. In section EJ we define a normalization 
function and extract its convergent factor, which gives us good representations of the ex- 
pectation values of arbitrary local operators in the ground state. We obtain a recursion 
relation for the normalization function which is useful to prove important properties of 
the observables. In section we evaluate the expectation values of local operators in the 
infinite- volume limit. We prove the cluster property of the domain wall ground state on 
the basis of the recursion relation. 

2 The Model 

2.1 Definition of the Model 

We define a one-dimensional lattice A with a "length" L by 

A:=jx = | y £ Z with \y\ < l| , (1) 

where L = 21 + 1 is an odd integer. A can be decomposed two sub lattices A G and A' 
defined by 

A D := {x G A | x is an integer} (2) 

and 

A' := {x G A | x is a half-odd-integer} . (3) 

We define a deformed flat-band Hubbard model on the lattice A. We denote creation 
and annihilation operators by cj. a and c XjCT respectively, for an electron at site x G A with 
spin a =t, I- They satisfy the standard anticommutation relations 

\_C x ,cri C l,r} = dx,y3<T,T> {Cx,ai Cy,r} = { C x,ai C l,r} = (A) 
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for x, y G A and a, r =f, |, where {A, B} = AB + BA. The electron number operator is 
defined by 



whose eigenvalue gives the number of electron with spin a at site x. In the present paper, 
we use the open boundary condition. This condition is given by defining c Xj(T = for 
| x | > L/2. We consider a Hubbard model with a Hamiltonian 



H — H^op + Hi nt . 
The interaction Hamiltonian is given by 

H int = Uy~] n x ^n X)i 



(6) 



(7) 



with a positive number U which represents an on-site Coulomb repulsion energy. The 
hopping Hamiltonian is defined by 



o-=T,i zGA' 



-p(«r)/4, 



•J>(<0/4, 



where p(a) — 1 if a =j and p(cr) = — 1 if cr =| with t > 0. We can represent the hopping 
Hamiltonian also in a standard form with a hopping matrix 



-'-'hop / v / ,, S.J/hE.CT 1 !/,' 



(9) 



where 



/..I i i _ i . r , 

|g| 2 + \q\ 2 it x = y G A 

A 2 if x = y G A' 

if x = y _|_ 1 and y £ A 

Ag pM/4 if x = y - | and y G A ' 

(g/g*)PW/4 if X = ?/ + 1 G A 

otherwise 



(10) 



The Hilbert space of the model is spanned by the following basis 



\xeA 



n 

Ky&B 



where A, B are subsets of A and $ vac is the unique normalized state with no electrons 
defined by c XjCr $ vac = for any x G A and a =], |. In the present work, we consider the 
Hilbert space 7i with the electron number fixed to L = |A G |, which is spanned by the 
basis states (JUJ with \A\ + \B\ = L. 1 

1 Throughout this paper, |„4| denotes the number of all elements in a set A. 
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t(kl'+lQl t(\q\? + \q\ 




t(\q\? + W\ 5) t(M' + kl ') 



Figure 2: Arrows and gray circles depict hopping amplitudes and on-site potentials respec- 
tively for electrons with up(a) and down(b) spins on a unit cell. For complex q = |g|e* 9 , 
the hopping amplitude opposite to the arrow is the complex conjugate of that along the 
arrow. 



2.2 Symmetries of the Model Hamiltonian 



Let us discuss the symmetries in our model. The local spin operator Sx and the total-spin 
operator are defined by 



^x^ n ^ ] C x,ctCP^)<t,t c x 



(12) 



and 



cU) _ V 9^ 



with j = 1, 2, 3, where V^' is the Pauli matrix 



1 

1 



% 



7?(3) 



1 
-1 



(13) 



(14) 



The first important symmetry is U(l) symmetry. The Hamiltonian commutes with 
S^l, and an eigenvalue of S^l is called a magnetization which takes ±1/2, ±3/2, 
±L/2. We define U(l) transformation by a following map U v which is generated by S. 
in a finite system: 

= e^ p(CT) 4,. and ^(c x , CT ) 
where is a real parameter. This map is linear 

U V (A + B)= U V (A) + U^B) 
and homomorphic 

U v (AB) = U v (A)U ip (B), (17) 

for arbitrary operators A and B. Each spin operator is transformed to the following form 
under the U(l) transformation: 



(3) 
tot 



(15) 



(16) 



Uy{s!»)=<xx<ps!p- 

U V (S®) = sin up + 



sin ip 

sin 90 b^' + cos</? 



c(2) 

C(2) 



(18) 
(19) 
(20) 
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Note that the U(l) symmetry is enhanced to the SU(2) in the case of q — 1, i.e. the 
Hamiltonian commutes with the all total-spin operators S^l (j = 1, 2, 3). 

The second important symmetry is a Z2 symmetry defined by a linear and homomor- 
phic map II: 

n(c a;jCT ) = c-z.a and n(4 >(T ) = cl^ (21) 

where a takes 1 if a =t an d t if cr =1. Note the Z2 transformation of the total spin 
operator n(S' tot ) = — 5* tot . An energy eigenstate with a magnetization M is transformed 
by II into another eigenstate with a total magnetization —M, which belongs to the same 
energy eigenvalue. 

The Hamiltonian of a system with a large volume is invariant under a translation 

T u (c x<J ) = c x + UjCT and T u (c XIJ ) = c x+u ^ a , (22) 

except several terms near the boundaries, where u is an integer. Thus, the translational 
symmetry appears after taking the infinite- volume limit. 



3 Ground State of the Model 

In this section, we construct ground states of the deformed flat-band Hubbard model. 
The construction method is based on the original method obtained by Mielke and Tasaki 
for the positive semi-definite Hamiltonian in the flat-band Hubbard model IE] . 

3.1 Localized Electron Operators 

We define localized electron operators a x a and d X)(7 by 

-<f (<T)/ V ! + A4 CT -<r p(CT)/ V , ifxeA 

if x e A' ' 




and 



2 ' 1 2 ' 



c x a and c x>a can be written in terms of a x a and d x y. 



and 



(23) 



d xa = { \° x ' a ifxGA D _ 
q~ p{ ° )/A c x i + Ac x , CT + q^ )/A c x+ , if x E A' 



ig p M/V , +ic4 CT + <r p(fT)/ V t ifxGAo 
4,a = <! x ~° A , (25) 

Aat „ if x e A' 



Xd X)(T if x G A 

-f Ka)/4 Ci, + ^ - <f . if x g A' (26) 

For an open boundary condition, we define c x<a = for \x\ > L/2 which means a xa — 
and d xa = for \x\ > L/2. The relations flU)), (|2U) , fl23) and (jUJ) are also valid 
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consistently under this condition. The inverse relations (J25|) and (J26|) implies that the 
Hilbert space Ti can be spanned also by the basis states 



®'a,b= n<r [u<i (27) 

where A and B are subsets of A with + \B\ = L. Note that a x a and d x>a satisfy the 
canonical anticommutation relations 

{4,a» d v,r} = S x ,y5 a , T , {4 )ff , aJ T } = dy, T } = 0. (28) 
The hopping Hamiltonian can be written in terms of d X)Cr , 

^hop = E <<A" (29) 

This form of the hopping Hamiltonian indicates the positive semi-definiteness i^hop > 0. 
The positive semi-definiteness of the interaction Hamiltonian H int > is also clear because 
n x ,a = c' x a c Xt a > 0. Thus, the total Hamiltonian is also positive semi-definite 

H = H hop + H int > 0. (30) 

First, we define all-spin- up and all-spin-down states by 

$ T= (n^Pvac and $ i= ( n<|) $ vac (31) 

\a;GA / \xeA / 

We easily verify H<&^ = = HQ i from the anticommutation relations of a x a and d X)Cr , 
and therefore these states are ground states of H. 

Next, we determine all other ground states. The conditions on a state $ to be a 
ground state H hop $ = and H irit $ = are equivalent to 

d x , a $ = (32) 

for all x G A' and 

CyjCy,^ = (33) 

for all y G A with a =f, [. 

We consider a restricted Hilbert space Hm spanned by the states with a magnetization 
M. Since the magnetization is a good quantum number, we find other ground states in 
each restricted sub-space TCm- Let $(M) be a ground state with a magnetization M. We 
expand $(M) into the following series 

$(M) = X>GA,S) ( ] (114,1 J $ vac, (34) 

.4,8 W.A / VyeB / 

where the sum is taken over all subsets A, B C A with \A\ = (L + 2M)/2 and |jB| = 
(L — 2M)/2, and ^(.A, £>) is a coefficient. The first condition (JH2*|) implies that ^(A., £>) 
does not vanish only for A PI A' = and SflA' = 0. The second condition (J3*3*j) for y G A G 

becomes 

c^A^M) = 0, (35) 



S 



due to the definition of d yjlJ . To satisfy this condition, tj)(A, B) must vanish for AC\B ^ 0. 
Thus the expansion (}3"4*|) can be written in the form: 

$(M) = II °U*w: ( 36 ) 

where Sm is a set of all possible spin configurations a = (a x ) x eA with a magnetization 
M: 



<Sm = \ ° — ( <J x)xeA 



^(-1)^=M, a,= U . (37) 



2 

XGA d 



In this representation, the second condition (J33|) for y G A' is equivalent to 

(4,t4+i,1 - ?4,|4+i,t) $ = (38) 

where we use inverse relations (I26|) and anticommutation relations (|28j) with x = ?/ — | G 
A G . We find that this condition is satisfied if and only if the coefficients satisfy 

0(a) = g [p( - ) - p(CT ^ l)]/2 0(a^ +1 ) (39) 

for x G A G with x < (L — 3)/2, where cr^+i is a spin configuration obtained by the 
exchange a x and c^+i in the original configuration a. This relation implies the uniqueness 
of the ground state with a fixed magnetization, since two arbitrary spin configurations 
can be related by successive exchanges of two nearest neighbour spins. Therefore the 
degeneracy of those ground states is exactly the same as that in the SU(2) symmetric 
model, as in the domain wall ground states in the XXZ model £Q. 
If we define S~ by 

The ground state with a magnetization M is obtained by acting this operator L — 2M 
times on the all-spin-up state 

$(M) = (^) L ^%. (41) 



3.2 Another Set of Ground States 

To explore the property of the ground state, we write it in a more explicit form as obtained 
by Gottstein and Werner in [T^]. We define the ground state ^(C) by 

*(0 = ( II «K0 J $vac = C(S-) n ^ (42) 

VzSAo / n=0 

where c4(C) := a' x * + (q x a x ,. We fix the order of the electron operators depending on 
their sites. In the XXZ model, the state which corresponds to this ground state ^(C) 
is a pure state in the infinite- volume limit, and therefore we can expect that this state 
is stable in the infinite- volume limit ^7]. Note that the state ^(C) defined here in the 
Hubbard model is not a product state unlike the corresponding state in the XXZ model. 
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3.3 Spin One-Point Functions 

A spin expectation value of a localized spin at site x in the ground state \&(C) is written 



as 



/can \- (T U) Ut (*(0, 4,^,^(0) \- (v^Ur MCW^j 

(S^)c= 2^ —o IKiTTTTl^ = n i^TTTJT, • (43) 



2 IWOII 2 
The following anticommutation relations with y G A 

{c XjCT ,aJ(C)} = A{4,T>4,t + = X Vx,Jx,y, 



mew 



for x G A„, and 



(44) 



(45) 



for x G A' are useful to calculate the expectation value, where rj x ^ = 1 and rj x ,i = Cq x - 
These anticommutation relations f)44j) and (|4"H|) yield the following equation 



c x ,*y(0 = sgD.[x]r] Xt „V x (C), 



(46) 



where ^(C) is defined by 



if x G A 



**(C) = < 



(47) 



The ordering manner of electron operators determines sgn[x] = ±1, which is irrelevant to 
expectation values of any local operators. Then, the expectation value of c x a c XjT for all 
x G A in the ground state ^(C) can be written in terms of ty x ((), 



mow 1 



,I*(C)II 2 

Also a density one-point function (n x )^ can be written in 



(4£ 



K) c = (i + iozT) 



..^>^CC) II 2 



l*(0ll : 



(49) 



where n x = n Xj | + n^j. Thus we obtain the representations of spin one-point function at 
site x in terms of a density one-point function (n x )c 



2 i + |Cg x | 2 ' 

(n x )((q x - {(q x )* 
2i 1 + |Cg x | 2 ' 

K) c i-ICgT 

2 l + |Cg1 2 ' 



> . 



(50) 
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Figure 3: The density one-point function (a) and the spin one-point function (b) are 
obtained by numerical calculation of rescaled normalization function for L = 101, q = 1.2, 
C = 1.2 -20 and A = 1.25. The white circles plot the values on A and the black dots plot 
the values on A'. The lines in (a) are asymptotic values of density one-point function 
shown in Eq. (JoTj) for A G and A' respectively. 



Here we show briefly the estimates of the density one-point function. We expect the 
density one-point function in the ground state ^(C) takes almost a constant value in A G 
or in A'. Indeed, we will prove this conjecture by evaluating the bounds of this function. 
Here, we show only the asymptotic forms of the function for large \x — z\. The details of 
evaluating the bounds are given in appendix IA.31 The asymptotic forms of the density 
one-point function are 

[ ,f [l + 0(p" 21 ^')] for xeA 

lim (n x ) c = { f " 4 2 v , (51) 

Ul-_£=J [l + 0(p- 2 '— I)] for xeA> 

where e = A 2 + |g|3 + \q\~z and p = min{r, |g|} for \q\ > 1 or p = r for \q\ = 1 with 
r = (e + V e 2 — 4)/2. In these estimates, the density one-point function is shown to be 
finite and strictly positive at each site for any finite parameters. Since the density one- 
point function is nearly equal to the asymptotic value, the behavior of the spin one-point 
functions is not controlled by the density one-point function besides the oscillation between 
integer sites and half-odd-integer sites. The spin one-point function (•S'i 3 ^ vanishes only 
at x = — logui | C|- As discussed for the domain wall ground state in the XXZ models 
El EI, the two domains are distinguished by the sign of the local order parameter 
(Sx )f. If the domain wall center is defined by the zeros of (S x /£, then the location of 
domain wall center is unique. The function \{n x )(^ — \ {Sx)c\ decays exponentially as the 
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function of the distance from the center. The domain wall width 1 / log \ q\ is defined by its 
decay length. In the large A limit in this model, the density one-point function is exactly 
equal to 1 for integer sites and for half-odd- integer sites. In this limit and for real q > 1, 
electrons are completely localized on integer sites, and the spin one-point functions are 
exactly the same as those obtained in the XXZ model defined on A D . 

For a complex q = \q\e l9 with a real number 8, one can see the spiral structure with 
a pitch angle 9. The vector (S x ) is rotated with the angle Ox around the third spin axis 
depending on the site x. Note that this spiral structure of the ground state does not 
exist in the XXZ model. Though the complex anisotropy parameter q = e %e is possible in 
the XXZ Hamiltonian, the spiral state is no longer a ground state and the corresponding 
model is described in the Tomonaga-Luttinger liquid without ferromagnetic order. The 
translational symmetry in the infinite-volume limit is broken by the domain wall or the 
spiral structure for a finite log |£|. Both symmetries U(l) and Z2 are broken spontaneously 
as well. We will discuss the broken symmetries in subsection 15.21 

3.4 Electron Two-Point Function 

An electron two-point function {c xa Cy jT )^ can also be written in terms of ^(C) an d ^y(C) 
by 

, t c , _ _ ( _ 1)y - x * n (tt*(C),tt y (0) (52) 

\ c x,a c y,r/C - ||\I>(C)|| 2 'lx,cr'ly,T \\^f(Q\\ 2 ' K ' 

We will see that the electron two-point function can be written in the form: 

( 4 £\ \v-A 

where K(x, y, L, () is a slowly varying finite function. This expression shows that the 
electron two-point function decays exponentially for large \x — y\ with the correlation 
length 1/ logr. 

3.5 Spin Two-Point Functions 

The spin two-point functions can be written in the form of 

(S^S?\= £ (^(^^(4 i<r ^^> f . (54) 
To obtain more explicit form, first, we rewrite the bracket into 

, + f v _ (Cz,o-Cy,o-'^(C); c x , T Cy , T '^ {(,)) 

\ C x,a C x,T C y,a' C y,T'/t — 11^(0 II 2 ' 

Using the anticommutator (JUj) and (fUjj). we can pick out the spin dependence as coeffi- 
cients 

c x , a Cy,T^(C) = sgn[x,y]r] x ^r] y ^ Xt y(C), (56) 
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where sgn[x, y] = ±1 and 
A 2 



n «uo 



A 



A 



n «-(o f -^ n 

w^x,y— i 



if x, y G A Q 
$ vac if x G A , y G A' 

$ vac if y G A OJ x G A' 



(57) 



MC) f - MO 1 



$ vac if x, y G A' 



Then the bracket in (|54jl is reduced to 



Note that the density two-point function is written 



ll*(C)ll 2 



(5? 



Kn,) ? = (l + |Cg1 2 )(l + |C^ 2 )^^. 



(59) 



We then obtained the representation of the spin two-point function in terms of the density 
two-point function as well as the spin one-point function written in terms of the density 
one-point function 



(60) 



cr.r.a' ,T 



+ \(q x \ 2 1 + IC^I 2 

If A is finite, we can rewrite spin two-point function from the spin one-point function (150(1 

(SSPSM)! = (Si\(S^) c {U ^\ . (61) 

K)cK)c 

In section EJ we will find that (n x )^ and (n x n y )^ exist in the limit |x| — »■ oo, |y| — ► oo and 
both. And we will prove 



namely, 



lim lim Un x n v )^ - (n x )^(n y )A = 0, 

\x— y\ — >oo L— >oo 



lim lim , = 1. 



(62) 



(63) 
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Thus, we conclude that the spin two-point function satisfies the requirement of the cluster 
property. We can also find bounds of the density two-point function. The function behaves 
asymptotically 



lim \{n x n y ) ( - (n a; ) c (n 2/ ) c | = < 0(p 2lv z| ) if \y - x\ > \y - z\ > \x - z\ (64) 



L— >oo 



\y 


— x\ 


> 


\x — z\ 


> 


\y - A 


\y 


— x\ 


> 


\y - A 


> 


\x — z\ 



0(jp 2 \v x \~j otherwise 



where p has the same value as that appeared in Eq.(foT|). We obtain these decay relations 
in appendix IA.21 The spin two-point function is almost a product of two spin one-point 
functions. This fact shows that the electron has a good localization property. In other 
words, \&(C) is almost a product state. In the large A limit, (n x n y )^ — (n x )^(n y )^ becomes 
identically, and then the spin two-point function becomes just a product of two spin 
one-point functions. This property is the same as domain wall ground state in the XXZ 
model. We believe that our model is related to the XXZ model if q G R. 

4 Normalization Function 

In this section, we introduce a normalization function which is useful to evaluate the 
density one-point function and two-point functions in the ground state ^(C)- 

4.1 Definition of Normalization Function 

We define a normalization function A(x, y\ () for x, y G A D with x < y by 

y 



A(x,y;C) 



I«I(C)K 



(65) 



Note that the norm of the ground state \I/(C) can be represented by ^(Oll 2 = ^4( — U h C); 
where I = (L — l)/2. The anticommutation relation for a x (Q and aJ,(C) for x, y G A is 

K(C), 4(0) = (1 + IC<f _l I 2 )e~ l k-i,, + e(l + |0f | 2 )4,, + (1 + \(q x+1 A 2 )e l ^ x+1 , y (66) 

where e = A 2 + \q\s + |(/|~^. Using this anticommutation relation, we can obtain the 
recursion relation for A(x, y \ C) in terms of A(x, y — l;() and A(x, y — 2; (), or A(x+1, y \ Q 
and A(x + 2,yX): 

A(x, y; z) = e(l + \(q y \ 2 )A(x, y - 1; C) - (1 + |C<f ^ | 2 ) 2 ^(x, y - 2; C), (67) 

or 

A(x, y; z) = 6(1 + \(q*\ 2 )A(x + 1, y; () - (1 + \ 2 f A(x + 2, y; ()■ (68) 

This forms show that we can obtain A(x, y;() recursively from the boundary values 
A(x, x; () and A(x, x + 1; (,), or A(y, y; () and A(y — 1, y\ (), which are easily calculated 

as 

A(x,x;C) =e(l + \(q x \ 2 ) 

A(x,x + 1; c)=6 2 (i + |C<f| 2 )(i + \(q x+1 \ 2 ) - (i + \(q x+1 A 2 ) 2 } ■ (69) 

A(y - 1, y; C)=6 2 (l + |Cg^T)(l + K<f | 2 ) - (1 + |C^^| 2 ) 2 
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The following formal boundary values are convenient 

A(x,x-l;() = A(x + l,x;C) = l, A(x,x-2;C) = A{x + 2,x;() = 0. (70) 

We can also represent A(x, y; () in terms of the values A{w + 5, y; () and A(x, w — 5;() 
for w G A D and 5=1,2 with x < w < y, 

A(x, y; () = -(! + |Cg-i \ 2 ) 2 A(x, w - 2; ()A(w + 1, y; C) 
+ e(l + |C^| 2 )A(x, w - 1; C)A( W + 1, y; C) 
- (1 + |0f^ | 2 ) 2 A(x, W - 1; ()A(w + 2, y; C), 



(71) 



where we can use the formally defined boundary values (j7()j) . Note that the relations 
(JfiTJl . (jHEJ) and (|7T|l are equivalent to each other with the formal boundary values (ffO|) . 
Hereafter, we consider that A(x, y; () is a function of x, y G Z and ( G C with y — x > —2, 
which is defined by the recursion relation (joTj) and formal boundary values ()70|) . 

If |g| = 1, then the coefficients of the recursion relation (joTj) are independent of y. 
This condition gives a great advantage for calculations. We define A(x,y) by 



A(x,y;C) = (l + \(\ 2 ) y - x+l A(x,y). 



(72) 



The recursion relation and boundary values for A(x,y) are obtained from those for the 
original function A(x, y; (), 



A{x, y) = eA{x, y - 1) - A{x, y - 2), 



(73) 



and 



A(x, x-2)=A(y + 2, y) = 0, A{x, x - 1) = A(y + 1, y) = 1. (74) 
By rewriting this recursion relation to 



A(x, y) — rA(x, y — 1) = r 1 A(x, y — 1) — rA(x, y — 2) 



(75) 



where r = (e + \/ e 2 — 4)/2, we can easily obtain the solution 

r 2/-a;+2 _ r ~(y-x+2) 



A(x,y) 



-i 



(76) 



Then we obtain a compact representation of the function 

, 1 1 _ r -2(y-x+2) 



A(x,y;C)= [r(l + |C| 2 )] ; 



(77) 



4.2 Rescaled Normalization Function 

In the case of |g| ^ 1, we parametrize \(\ by a real number z as |C| 
rescaled normalization function B by 



|g| z and define a 



\2(w-z) 



) 



B(x,y,z). 



(7t 
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The corresponding recursion relation for the rescaled function obtained from those for the 
original function (|57j) is written 

1 + r 2 1 — f(v — z — -) 
B(x, y, z) = —^B{x, y - 1, z) ±L B {x, y-2,z) (79) 

where / is defined by 



\2 



/(«) = 1 i i -• (80) 

+ |g| tt -3)(|g|- u -3 + [g|«+3) 

This /(u) is an even function of a real number it. At u — 0, /(w) takes the maximal value 

/(0) = f | g |!~| g f^ , (81) 
\|g|a + |g| 2y 

and for large \u\, f(u) is proportional to |g|~ 2 l"l asymptotically. The boundary values are 

B(x,x-l,z)=B(x + l,x,z)=l, B(x,x-2,z) = B(x + 2,x,z) = 0. (82) 
The relations obtained from (|68|) and ()71|) are 

1 + r 2 1 - fix - z+ -) 
B(x, y, z) = ~^B(x + 1, y, z) ^5(x + 2, y, z) (83) 

and 

B(x, y, z) = ^— ^ w - 2, z)B(w + 1, y, z) 

1 + r 2 

-B(x, w — 1, z)B(w + 1, y, z) 



/(w + | - z) 

2 



B(z,iu - l,z)-B(iy + 2,y,z), (84) 



which are equivalent to the recursion relation (|79|) . 

Let us discuss properties of the rescaled function. The function B(x,y,z) is non- 
negative and monotonically increasing in y — x, i.e. 

B(x,y,z)>B(x,y-l,z), B(x, y, z) > B(x + 1, y, z). (85) 

To obtain these relations, we rewrite the recursion relation (|79|) 

p/ \ r/ in B (x,y- l,z) -B(x,y-2,z) f(y-j-z) 

B{x, y, z) - B(x, y-l,z) = = + j 5(x, y - 2, z). 

(86) 

This representation shows B(x, y, z) > B(x, y — 1, z) if -B(x, y — 1, z) > -B(x, y — 2, z) and 
B(x,y — 2, z) > 0. From the formal definition (j%2~|) . we know j?(x. .x — 1. z) > B(x,x — 2) = 
, and therefore B(x,y,z) > B(x,y — l,z) > 0. Also B(x,y,z) > B(x + l,y,z) is 
shown almost in the same way. The function B(x, y, z) is invariant under replacing |g| by 
|g| _1 since the coefficients in recursion relation (fTTJj) and the formal boundary values are 
invariant under the replacement of |g| by |g| _1 . 

Though we cannot obtain a compact form of the solution B(x,y,z), we show that 
B(x, y, z) exists in each limit y — > oo, x — > — oo or both respectively in the next subsection. 
And we can obtain bounds of the ratio of rescaled normalization functions at two different 
sites. This bounds are useful to estimate the expectation values in the ground state ^(C)- 
We will obtain these bounds in appendix I A. 11 
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4.3 Convergence of the Rescaled Function 

Let us prove that B(x, y, z) converges in each limit y — > oo, x — > — oo or both. We use the 
relation (}84|) for these estimates. The existence of _B(— oo, oo, z) is shown by this relation, 
if we show the existence of B(x, oo, z) and B(—oo,y,z). We then concentrate to show 
the convergence in the limit y — > oo. The convergence in another limit x — > — oo can be 
shown in the same way. 

We introduce new functions C(x,y,z) and D(x,y) where x, y G Z and 2 e M with 
y — x < —2. These functions are defined by following recursion relations: 

L>0, y) = ^p^D(x, y-l)- ^-D(x, y-2), (87) 



and 



1 + r 2 1 

V, z) = — 5— C(x, y-l,z)- -zC(x, y-2,z) 



9 "D3 I 9 



T ^—^D{x,y-2). V oo; 

with initial values 

C(x,x-2,z) =D(z,x-2) =0, C{x,x-l,z)=D{x,x-l) = l, (89) 

where (7 is a real parameter. Note that both functions C(x, y, z) and D(x, y) are nonnega- 
tive and monotonically increasing for y — x as well as the rescaled normalization function 
B(x,y,z). Both recursion relations for C(x,y,z) and D(x,y) are much easier than that 
for B(x,y, z). One can easily obtain a compact representation of the solution D(x,y) 

jjy-x+2 _ py-x+2 

D(x,y)= K+ R+ _ R (90) 

where R± are defined by 

R + = 1 + 7(0) and R_ = r~ 2 - 7(0) (91) 

with 

v /6 2 -4 + %-y^4 
7(5) = ^ • (92) 

Note that j(g) takes a range < j(g) < r~ 2 for < g < 1. 
If g satisfies 

<7 > max /(ty + i — z), (93) 

x<w<y— 1 

these functions satisfy the following inequalities, 

B{x,y,z) <C{x,y,z) < D(x,y), (94) 

and 

2/, z ) - B (x, y-l,z)< C(x, y, z) - C(x, y-l,z), (95) 
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for any integers x, y with x < y and a real number z. To show the second inequality in 
flUIJ), we calculate D(x, y) — C(x, y, z) using the recursion relation for each 

D(x, y) - C(x, y, z) - [D(x, y - 1) - C(x, y - 1, 2)] 
= i{ D(x, y-l)- C(x, y-l,z)- [D(x, y-2)- C(x, y-2, z)} } 

+ fiV ~J ~ Z) [D(x, y-2)- C(x, y - 2, z)] + 9 - fiV ~ I ~ " } D(x, y-2). (96) 
And note 

a — f( x ±I_ 2 ) 

D(x, x) - C(x, x,z) = and D(x, x + 1) - C(x, z + 1, z) = - — ^ 2 2 > 0. (97) 

Eqs. (JOOJ) , (j^Tjl and the condition (jHSJ) imply that D(x, y) —C(x, y, z) is non-negative and 
monotonically increasing function of y — x > because D(x,y) > 0. Thus, we obtain 
D(x,y) > C(x,y,z). The remaining two inequalities in (}9"4"j) and (|95|) can be shown in 
the same way. 

Now we prove the convergence of C(x, y, z) for g = /(0) in the limit of y — > 00. The 
convergence of C(x, z) and the first inequality in (|94|) give the convergence of -B(a;, y, z). 
We rewrite the recursion relation for C(x, y, z) 

n t \ n < -\ \ C(x,y-l,z)-C(x,y-2,z) 
C{x, y, z) - C(x, y-l,z) = — 2 



f{y-\-z)Rr-RT x 

+ ? R+-R. ■ ( ® 



The final term can be evaluated by 
f{y-\-z) R y - X -R y ~ x 

r 2 R+ - R_ 

R_ \" 



1 










1 

2 — 


q\ 


1 

2 


? 




R -x _ R v_-y R v + 




q 


-22-1 _|_ 




-Ay+2z+l 


+ 


\q 


-2y-l + 


q 


R + - R_ 



< 



R+Y \q\ 2z+1 (\qfi 



1 1 —i 
\q\ 



\q\ 2 J r 2 R\ R + - R^ 
From the recursion (|98|l. we have 

nt \ nt 1 \ ^ C(x, y — 1, z) — C(x, y — 2, z) Q n T s t \ tnn\ 
C{x, y, z) - C(x, y-l 7 z)< ~ 2 + fi v K x {x, z) (99) 

2 



where (3 = R + /\q\ and 



22+1 " 



* (\q\ \ - 



^^^W' < 100 > 



Here, we show (3 < 1 for \q\ > 1 and g < /(0). We define a function 



h ^ g > = , /~2 — 7 ■■ 

e + v e — 4 
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which is monotonically decreasing in e > 2 and monotonically increasing in g > 0. Then, 
we represent (3 as 



\q\ 

Since e = A 2 + \q\z + \q\~^ >2,g< /(0) and \q\ > 1, we obtain 
P = -Lh(e ig )<^h(2,f(p)) = 



3/1,1 I, 

2 + m 



TT<1- 



:ior 



(102) 



Using (|99|) recursively, we obtain 

C{x, y, z) - C{x, y-l,z) < - A + —^-K^x, z) + P y K 1 (x, z) 



where 



and 



1 (r 2 B)~( y ~ x+ V — 1 

=K 2 (z - x)r-^ y - x) + K 3 (z - x)(3 y - x 



K 2 {z - x) 



KJz - x) 



\q 


]2(z-x)+3 n q 


i 

2 — 


Q 




\q\ 


2 - r 2 R + R + - 


R- 



R + 


q 


2(z- 


x)+l (| ? 


i 

2 — 


q 




r 2 R. 




q 


2 R + - 


R. 



(103) 
(104) 

(105) 



-^2,3 (2 — 2) are finite for fixed x and z. Then we have found that the difference between 
C(x, y, z) and C(x,y — 1, z) decays exponentially as y becomes large if x and z are fixed, 
and hence C(x,y,z) converges for y — > 00. This fact means that the limit B(x,oo,z) 
exists. As discussed before, we easily find the convergence in the limit of x — > — 00 and 
in the both y — > 00 and x — > —00. 



5 The Infinite- Volume Limit 

In this section, we consider the infinite- volume limit. First, we define a ground state 
uj^ in the infinite- volume system on the basis of the ground state ^(C) i n the finite- 
volume system. We show that any local operator cannot deform the ground state uj^ to 
a state with lower energy than that of Uq. Next, we remark the symmetries of the model 
Hamiltonian in the infinite-volume limit and discuss their breakdown. We show also the 
cluster property of the ground state and excited states above uj^. For the general proof, 
we use the representations in terms of the rescaled normalization function which is defined 
only for \q\ 7^ 1 in this paper. For \q\ = 1, however, the proof is much easier. 

5.1 Ground State in the Infinite- Volume System 

A state co>£ is a linear functional over all local operators for the infinite-volume system 
defined by 

u ( (0 A ) := hm(0 A ) c = hm , (106) 
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for an arbitrary local operator A written in electron operators on a finite subset A 
of the lattice A. We define also other states and uj^ for the infinite-volume system 
corresponding to the finite- volume ground states $| and $^ with the definition (|31|) . 
This state is the ground state in the infinite-volume system, if the inequality 

lim u c (O A [H Ll ,O A })>0, (107) 

is valid for any local operator A , where Hl 1 is a local Hamiltonian defined on a finite 
chain with length L±. This inequality implies that any local operator cannot decrease the 
energy, and hence we can regard the state as a ground state in the infinite-volume 
limit. In the Hamiltonian @, there exists L 2 G Z for the set A such that 

lim [0 A ,H Ll ] = [0 A ,H L2 }. (108) 

Li— >oo 

Since Hl 2 is positive semi-definite and Hl 2 ^(C) = f° r an y finite L 2 , 



mew 



> (109) 



is trivial. Thus, satisfies the inequality (|107|) . We can easily verify that the states 
and u>i are the ground states of the infinite-volume system. 

5.2 Symmetries and Their Breakdown 

Here, we consider the U(l), the Z2 and the translational symmetries of the Hamiltonian in 
the infinite- volume system. These symmetries are defined for a finite system by Eqs. fTo^ . 
(|21|) and ([22)1 . Let Q be one symmetry group of them. An arbitrary transformation X G Q 
can act on an arbitrary local operator A written in electron operators on a finite subset 
A of the lattice A 

A ^X(0 A ). 

The Hamiltonian of the infinite-volume system has a symmetry Q, if a local Hamiltonian 
Hl 1 defined on a finite chain with a length L\ has the following transformation property 

lim X([0 A ,H Ll }) = lim [X(0 A ),H Ll ] (110) 

Li— >oo Li^oo 

for any X G Q and for any local operator A . Thus, our model has U(l), Z 2 and 
translational symmetries defined by the mappings (|15|). (|2ip and ()22|) in the infinite- 
volume limit. 

The spin one-point functions in the infinite-volume limit are shown 



-c (^ 2) ) 
-c (^ 3) ) 



2 1 + \q x ~ z \ 2 ' 
U({n x ) q x ~ z - (q x - z )* 
2i 1 + \q 

uJ ( (n x ) 1 - \q 
2 l + |<f- z | 2 ' 



x-z\2 



x-z\2 



iir 
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All the symmetries of the model are spontaneously broken in the ground state u^, since 

vdU v (SP)) + u; c (U(S^)) + Wf (S?)), w c (T u (S< s >)) ? Wc (Sf ). (112) 

On the other hand, the spin one-point functions in the all-spin-up and all-spin-down 
ground states are 

^(SP) = = u,(S x \ ^ 3) ) = l^(n x ), (113) 

^(S«) = = ^(Sf), ^(Sf ) = -\^{n x ). (114) 

Therefore, the ground states and uo^ have the U(l) and the translational symmetries. 
The broken symmetry is only Z 2 in and u^. 

5.3 The Cluster Property 

First, we give a precise definition of the cluster property. Let A and B be finite subsets 
of A. And we denote a translation of A by A + x for x G A: 

^ + x:={x + y|yG A}. (115) 

We say that a state has the cluster property, if the truncated correlation function of 
arbitrary observables A and B on arbitrary finite subsets A, B C A satisfies 

lim [u{0 A B±x ) - u(O a )u(O b±x )\ = (116) 

x— >oo 

and 

lim KCU ±:r O s ) - uj(O a±x )u(O b )} = 0. (117) 

Here we prove the cluster property of the ground state u^. Let AjB is a basis of local 
operators: 



o A ,B( a A, <tb) = ( n ^ ) ( n d v^y 



:ii8i 



where a A = (cr x ) xeA with a x =|, |. The relation given in (j^Bj) means that any local 
operator can be written in the form 

Yl Yl &(^>°* )(>A.K.(^A : ^t<:) (H9) 

where (j)j(cr Aj ,crB ) G C is a coefficient. Hereafter we do not write the spin dependence 
explicitly for simplicity. Note that the summation in this representation ()119j) consists of 
finite terms. The cluster property is shown if we find 



and 



lim 

x— >oo 



lim 

x— >oo 



^d°ABOe+x,v+x) - ^d°A,B)^d°c+x,v+x) =0, (120) 



WcSP A-x,B-xOc,T>) — ^((O A~x,B-x)^C,{Oc,v) 



o, ri2i: 



for any finite subsets A, B, C and T> of A. The following two lemmas are useful to prove 
these equalities: (THZUJ) and (fT2"Tj) . 
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Lemma 1 Let xj G A (j = 1,2, ••• , |v4|) and y^ G B (k = 1, 2, • • • , |£>|) be ele- 
ments labeled with the monotonically increasing property Xj < Xj + i and < yu+i- If 
{Oa^^A: °'b))c does no ^ van ish, then A and £> satisfy the following conditions: A,Bc A Q , 
|.A| = \B\, Xj < yj =>- yj < Xj + i and yj < Xj =>- Xj < yj + i for any j = 1, ■ ■ ■ ,\A\ — 1. In 
this case, (Oa,b)( can be written in the form of 

in (~ ~ \\ B{-l,u-l,z)B{v + l,l,z) 

{Oa,b{°a, vbMc = B(v + 1 I z) T AB\ a A, vb), (122) 

where u = min{xi,|/i} and v = max{xui, ym\\- Ta,b and its translation Ta± w ,b±w are 
finite even in the limit of \w\ — > oo. 

Lemma 2 We consider four finite subsets A, B, C and T> of A c , and define lattice sites 
x = min^4 U B, y = max A U B, u = minC U T> and v = maxC U T> with y < u. If 
{Oa,bOc,v)c does not vanish, we have to consider only the following two cases. If either 
{Oa,b)( or (Oc,v)c does not vanish, then it can be written in the form: 

ln n v B(-l, x - 1, z)B(y + 1,^-1, z)B(v + 1, 1 z) 

(yA,BUc,T>)c = j ^ J-a,bJ-c,v- (i-^o) 

The finite coefficient Ta,b is identical to Ta,b(&a> &b) i n Lemma 1. We have abbreviated 
the spin dependence. If both {Oa,b)c and (Oc,v)c vanish, the following function decays 
exponentially 

in O \ B(-l,x-l,z)B(v + l,l,z) ( e*lY~ y ~ 

(Oa,bOc,v)c B(-l,l,z) [~ ) (124) 

Ta,b,c,v and its translation Ta-w',b-w',c+w,t>+w is finite even in the limit \w\ — > oo, \w'\ — > 
oo and both. 

We can easily prove these lemmas in practical calculations of expectation values in the 
ground state ^(C) as for the density two-point function or the electron two-point function. 
They will appear in appendices I A. 41 and IA.5I 

Proof of the Cluster Property Now we prove the cluster property of the ground 
state ujq and excited states above ui^. In the case of {Oa,b)c = 0^ lemma 2 implies 
{Oa,bOc,v)( —> in the limit of (max A U B — minC U D) — > 00. Next, we consider the 
case of {Oa,b)( 7^ and {Oc,v)c 7^ 0- The lemmas 1 and 2 ensure that the truncated 
correlation function can be written in the form of 

(Oa,bOc,v)<;—(Oa,b)c(Oc,v}c 

B(-l,x-l,z)B(v + l,l,z) 

x [B(y + 1, u - 1, z)B(-l, I, z) - B(-l, u - 1, z)B(y + 1, 1, z)\ . (125) 

We already know that T c+WiT >+ w and lim^oo B{v + w + 1, 1, z)/B(—l, I, z) are finite even 
for w — ► 00 and 

lim lim (B{y + 1, u - 1, z)B(-l, I, z) - B(-l, u-1, z)B(y + 1, 1, z)) = 0. 

u — >oo I — >oo \ / 
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This relation enables us to show the requirements (|12U|) and (|121|) . Thus, the cluster 
property of the ground state uj^ has been proved. We can prove the cluster property 
of an arbitrary excited state above the ground state basically in the same approach, 
since any excited state can be obtained by acting some local operators written in a finite 
summation of the basis (jll8j) on the ground state uuq. In appendix IA.2| we will find the 
following bound for a sufficiently large u 

lim (B(y + 1, u - 1, z)B(-l } /, z) - B(-l, u - 1, z)B(y + 1, I, z)) < Cp~ 2(u " w) , (126) 

£— >oo V / 

with a positive number C independent of u for w = max{y,z}, y = max A U B and 
u = minC U T>. Here, we define p = min{|g|,r} for |g| > 1, or p = r for |g| = 1. This 
bound and the representation (J 124)) mean that all truncated correlation functions decay 
exponentially. And their decay lengths are less than (21ogp) _1 . 

The proof of the cluster property or other evaluations of expectation values for the 
ground states and are much easier than those for the domain wall ground state 
because the recursion relations of the normalization functions for $| and <3>^ are exactly 
solvable. 



Appendix 

Here, we show the estimates of the rescaled normalization function to evaluate the density 
one-point function and several correlation functions. In the case of \q\ = 1, explicit and 
compact representations of expectation values are easily obtained by the exact formula 
of normalization function (|77j) . Also in the case of |g| ^ 1, the dominant part of the 
two-point functions is the same as those for \q\ = 1. Here, we show the evaluation only 
for \q\ > 1 without loss of generality. 

A.l Estimates of the Rescaled Normalization Function 

Here, we show useful inequalities 

l-G + (x,y,z) < — — — r < 1, 1 -G-(x,y,z) < — < 1, A.l 

B{x,y + l,z) B{x-l,y,z) 

for x < y. The positivity 1 — G±(x,y, z) > is guaranteed by the definition 

G±(x, y, z) = min | ^ - , G±(x, y,z)\ , (A.2) 
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where the functions G±(x,y,z) are defined by 

r -2(y-x+l) 



G + (x,y,z) 



1 + r s 



, , r -2(j/-x+l) _ | fl |-2(i/-x+l) 
, ,_ 2 (x-*-§)_' 



14 



-2(z- 



+ (kl 2 - l?l ~ 2 ) 2 x < 



— r 



-2(3/-a;+l) 





g 2 _ r 2 




-2(j/-x+l) r -2( 2 /- a: +l) 




g 2 r 2 _ I 


|g|- 2( - 

o i In 


-x-A) ^1-2(2,-2+1) 



if a; > z — \ 



Xy<z-k 



(A.3) 



+ |g| r -2(»-' + l) 



r 2 |g| 2 — 1 |g| 2 — r 2 



otherwise 



and 

G-(x,y,z) 



-2{y-x+l) 



1 + r 2 



k/l 



-2(»-»-f) l ~ |g|- 2 fa"" +1 )r- 2 fa-" +1 ) 
\q\ 2 r 2 — 1 



-2(^-2/- 



^ r -2(y-x+l) _ | ? |-2(j/-x+l) 



+ (l<?l 5 " l<?l~ 5 ) 2 x < 



-2(j/-x+l) 



|g| 2 — r 2 

| 5 |_ 2 (y-*-i) | 9 |_2(*-a!+i) 



if x > z 



if 2/ < z + i 



(A.4) 



r 2|g|2 



r 2 — |g| 2 



m \r 2 \q\ 2 - 1 \q\ 2 - r 2 



otherwise 



where z' is an integer defined by z' G Z D [z — h, z + |). These inequalities give the 
estimates of the expectation values. Note that in a limit A — ► oo, r goes to oo, and 
therefore B(x,y — 1, z)/B(x,y, z) — > 1 for arbitrary x, y and 2. 

Now we prove the inequalities (jA.ljl . The upper bounds are obvious, because the 
function B(x, y, z) is monotonically increasing for y — x. The remaining task is to obtain 
the lower bounds. 

Here, we treat only B(x,y, z)j B(x,y + l,z). The lower bound for B(x,y,z)/B(x — 
1, y, z) can be obtained in almost the same approach. First, we show B(x, y, z)/B(x, y + 
1, z) > 1 — (r 2 + I) -1 - Using the relation (fTHJ) . we obtain 



B{x,y,z) 
B(x,y + l,z) 



1 + r 2 l-f(y+\-z)B(x,y-l,z) 
r 2 r 2 B(x,y,z) 



(A.5) 



Note B(x, y — 1, z)/B(x, y, z) > and 1 — f(y + \ — z) > 0. Then we obtain 

B(x,y,z) r 2 1 1 



£(:e,j/ + 1,z) 1 + r 2 



(A.6) 
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To obtain another lower bound, we rewrite the recursion relation 

I fly -|_ I _ z \ 

B(x, y + l,z)- B{x, y,z) = — [B(x, y, z) - B{x, y-l, z)\ H § B(x, y-l,z). 

(A.7) 

And using this relation recursively, we obtain 
B(x, y + l,z)- B(x, y, z) =— [B(x, y - 1, z) - B(x, y - 2, z)\ 



/• 1 



f{y-\- z ) ul n , + I c 

— - 1 y - 2, z) + 1 B(x, y-l,z) 



j=0 

by using B(x, x — 2, z) = 0. Then we obtain 



■2(y-x+2) + 2_. ^2(j+l) -1,^), (A.{ 



£(z,y,*) _ 1 1 1 g/(y + f-j-z)5(x,y-j-l : 



B(x,y + l,z) r 2(y-x+2) B(x,y + l,z) r 2 ^' +1 ) B(x,y + l,z) 

>1-- V 7 l ^ 2 1 fi (A 9) 

1 + r 2 r 2(i+l) 

by using 

1 + r 2 , B(x,y — j — l,z) 

B(x,y + l,z)>^- and „ , <1 A.10 

for y > x. Note the inequality 

f(u) < (|g|* - |g|~5) 2 |g|~ 2|tt| . (A.ll) 

We can evaluate the summation in the final line of (|A.9J) . There are three cases. First 
case is x > z — \. Since y + | — j — z > for < j < y — x, we obtain 

^ f(y + 1 - j - z) < (\ q \h - \ q \^y ^ 

j=0 j=0 

, , , r -2{y-x+l) _ \ \-2{y-x+l) 

= - ir^kr 2 ^- • (A12) 

For the second case y < z — - we obtain 

IJ-a: 



Z^ r 2(j+l) — r 2 Z^ r 2j 

i=o j=o 



i i , 1 - \ n \-2(y-x+l) r -2(y-x+l) 

(l«l* - kr^^l-^-'-i' 1 191 _ f t , (A.13) 
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because y + | — j — z < for < j < y — x. For the third case x < z — | < y, we 
introduce an integer n such that < y — n — z -\- ^ and y — {n + 1) — z + \ < 0. Then the 
summation can be evaluated by 

V-* *t.. , 1 l-l-i^r « | g |_2(»-j-»+^) ^ | 9 |2(j,_j_ 2+ I)- 

. j=0 j=n+l 



^ /(y + s-i-*) < (l^-kH) 



j=0 



r 2(j+l) 



(ki 5 - i?r 



r -2 fa -, +1 ) k|- 2( — g 

r 2 |g| 2 - 1 
+ |g|r- 2 ^' +1 > 



+ 



r 2i 

| ? j-2(y-z+i) 



|ff| S 



|g|2(«'-*) |g| 2 («-z') 



r 2 |g| 2 — 1 \q\ 2 — r 2 



(A. 14) 



where z' := y — n E Zn[z — \,z + \). Then we obtain another lower bound 



B(x,y,z) 
B(x,y + l,z) 



> 1 



r -2(y~x+l) 

1 + r 2 
\q\ 



(\q\2 - \q\ 2) 2 x < 



2(x- 




«) _ \g\-2{y-x+l) 






\q 2 — r 2 


2(2- 


y-l) 1 - 


. 2 (j / - a! +l) r -2(y-x+l) 






g 2 r 2 _ ]_ 


-2(1/ 




*-3) | g |-2(y-^+i) 



if x > z — k 



if y < z - \ 



(A.15) 



+ \q\r 



-2(y-z'+l) 



r 2 \q\ 2 — 1 r 2 — |g| 2 

/ [gp'- *) { |g| 2 ^') 



otherwise 



\r 2 |g| 2 — 1 \q\ 2 — r 2 

A. 2 Estimates of the Truncated Correlation Functions 

Here, we prove the inequality (|126|) . That inequality means that all truncated correlation 
functions decay exponentially. First, we show the simple form of bounds of derivation 
of two rescaled normalization function. For any x and z, there is a positive number Co 
independent of y such that for sufficient large y 



B(x,y + l,z) 



(A.16) 



with w = max{z,x}. Here, we define p = min{|g|,r} for \q\ > 1, or p = r for \q\ = 1. We 
can obtain this bound easily from the bounds ()A.1J1 . And following bound is trivial 

B(x,y + l,z) 



< 1 + C lP 



-2{y+l-w) 



B(x,y,z) 

Here, the finite positive number C\ is also independent of y. 
Let us prove the inequality (J126J) . 

B(-l,l,z)B(x + l,y-l,z)- B(-l,y- 1, 1, z)B{x + 1, 1, z) 

B(-l,y-l,l,z)B(x + l,l,z) 



(A.17) 



B(-l,l,z)B(x + l,y-l,z) 



B(-l,l,z)B(x + l,y-l,z) _ 



(A.18) 
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We concentrate to obtain bounds of 

B(-l,y -l,l,z)B(x + l,l,z) 
B(-l,l,z)B(x + l,y-l,z) 



(A.19) 



since we have already known that the pre-factor B(—l, I, z)B{x + 1, y — 1, z) is finite even 
in the limit of / — > oo and y — > oo. First we estimate an upper bound. Note that B(x, y, z) 
is a monotonically increasing function of y — x. We can evaluate 



B(-l,y-l,l,z) B(x + l,l,z) B(x + l,l,z) 



B(-l,l,z) B{x + l,y-l,z) B(x + l,y-l,z) 
And using Eq. (|A.17|) . we obtain 

B(x + l,l,z) B(x + l,l,z) B(x+l,l-l,z) B(x + l,y,z) 



(A.20) 



B(x + l,y-l,z) B(x + l,l-l,z) B(x + 1,1-2, z) B(x + l,y-l,z) 

i 

<Y[(1 + C 1P - 2{V - W) ). (A.21) 

v=y 

To estimate this, we take the logarithm. Then we can evaluate 



log 



2(v— w)\ 



Lv=y 



= ^log(l + C lP " 2 ^) 

v=y 

I 



(A.22) 



v=y 



for some positive numbers C 2 and C3 independent of y. Here, we have used the mean-value 
theorem to obtain the second line. Then we obtain an upper bound 



< gCiCp-^-p-'O-)) < l+c ^ v -2 { y- W )_ v -2 { l- w)) (A>23) 



B(-l,y- l,l,z)B(x + 1,1, z) 
B(-l,l,z)B{x + l,y-l,z) 

for a positive number C4. We also can evaluate a lower bound 

B(-l,y - 1,/, 2)5(0; + 1,/, ^ 



B(-Z,i,z) J B(x + l,y- l,z) 



> 1 - c h {v~ 2( - v ~ z) -p- 2{l - z) ), 



(A.24) 



with a positive number C5 independent of y. Then we have obtained the bound for a 
sufficiently large y 



lim [B(-l, I, z)B{x + l,y-l, z)- B{-1, y-l, l, z)B(x + 1, /, 



l— >oo 



< Cp- 2[y - w) . 

(A.25) 



A. 3 Estimates of the Density One-Point Function 

Now we can estimate the density one-point function. We first estimate the factor 

B(-l,x- l,z)B(x + l,l,z) 



B{-l,l,z) 



(A.26) 
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which appears in the representations of the density one-point function both for integer 
sites flA.30|) and half-odd-integer sites (jA.31|) . Using the relation (jH3), we obtain 



B(-l,x- l,z)B(x + l,l,z) 
B(-l,l,z) 



1 + r 2 _ 1- f(x-±-z)B(-l,x-2,z) 
r 2 r 2 B(—l, x — 1, z) 

1- f(x + ±-z)B(x + 2,l,z) 
r 2 B(x + l,l,z) 

This form and the inequality (|A.1J) enable us to estimate the upper bound 
B(-l, x- l,z)B(x + l,l,z) 



by using 



and 



B(-l,l,z) 
B(-l,x-2,z) 



< 



r 2 -l ^f^V 



B(x + 2,l jZ ) < 1 



B(-l,x-l,z) ~ ' B(x + l,l,z) 



1 - f(x ± \ - z) < 1. 
For the lower bound, we obtain the inequality 

B(-l,x-l,z)B(x + l,l,z) fl + r 2 I- f(x-\- z] 



B(-l,l,z) 



> 



r 1 
: V 4 2 ^41 + G 1 (x,/,^)' 



[l-G + (-Z,z-2,*)] 



[l-G-{x + 2,l,z)] 



(A.27) 



by using 

B{-l,x- 2,z) 
B(-l,x- l,z) 

with 

Gi(x, I, z) 



>l-G+{-l,x-2,z), B B ^ X \\YI - 1 ~ G - {X + 2 ' 1 ' Z) ' (A - 28) 



f(x-l-z)[l-G+(-l,x-2,z)} 



r 2 — 1 



+ \ -z)[l -G-(x + 2,l,z)] +G + (-l,x -2,z) + G_(x + 2,l,z) j. (A.29) 

The representation of the density one-point function in terms of the rescaled normal- 
ization function is 

A 2 B(-l,x -l,z)B(x + l,l,z) 



B{-l,l,z) 



(A.30) 
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for x G A G or 



191*"*+ M 



-x+z 



B(-l,x -^,z)B(x + l,l,z) 



\q\ s 



+ \q\ 



-X+Z + 



B(-l,l,z) 



\ q \x-z + \q\-*+z B(x + ll,z) 
r \q\*-*+l + \q\-*+*-l B(x + §,Z, z) 
| g |*-* + x-l z)B(x + §, Z, z) 



-x+z- 



B(-Z,Z,z) 



1 1 |g|*-' + | g |-*+* £(-Z,x-f,z) 
r | g |*-*4 + |g|~*+*+§ B(-Z, x-\,z) 



for x G A'. 

From the representation ()A.30j) . we find the bounds 



A 5 



1 



Ve234l + Gi(a;,Z,z) 



X 2 



(A.31) 



(A.32) 



Z, z) consists of terms which decay exponentially as Z — Z — \z\ or \x — z\ become 
large. And for sufficiently large \x — z\, lirn^oo Gi(x, Z, z) is proportional to p~ 2 \ x ~ z \ 



asymptotically: 



lim Gi(x,Z,z) = 0(p" 2|:c - z| ) 



(A.33) 



where p = min{|g|, r}. 

Next, we estimate the density one-point function for a site in A'. For upper bound of 
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the density one-point function, we evaluate 



\q\ x -*+\q{ 



|g|*-*-3 + \q\- x+z +l I r |g|*-*+3 + |g| 



l-i l 'l" +|g| ^'., [l-G-(» + iU)] 



+ 



\q\ x - z + \q\ 



-x+z 



\q\ x ~ z +i + \q\- x+z ~\ 



i- 1 - + l t i- C+ (-u-i»)l 



_|_ 1^1— x+2 7-|g| x ' — 2 +5 -)- r|g| _:E+z ~5 — [g^ -2 — |g| — x+z 



r(|g| z 



z+± 



■x+z- 



+ 



z +|g| x+2 r|g| a 2 2 + r|g| x+z +2 — |g| : 



-x+z 



1 |g| x_li + |g| _x+a 



■x+z- 



T G-{x + ±l,z) 



T G+(-/,x 



2 • 



\q\*- z + \q\ 



-x+z 



(\q\x-z-l + \q\~ x + z +k)(\q\*- z H + \q\- x + z ~l) 



Ve 2 ^! 

x I |?r*(kl* + I?! -4 - 2r" x ) + \q\- x+z (\q\ l > + \q\~i - 2r~ l ) \ + G?(x, I, z) 



1 



Ve 2 ^! 
1 



(v 7 ? 3 ! - A 2 )(|g|^ + \q\- x+z ) 2 



(\q\ x - z —2 + \q\- x+z+ i)(\qf 
[1 + Gt(x,l,z)i 



+ \q\- 



1 \ "I" (-^J ^1 



(A.34) 



where we have used \q\* + ]g| 2 = e — A 2 , and we determine a new function G\, which 
includes the terms which are proportional to |g|~' z ~ 2 '. And we obtain the lower bound 



(n x )c > 



1 



+ 

1 



| ? |*-* + \ q \-x+z 



1 |gr _2 + |g|- x+2 



1 



14 



r | g |x- 2 +| + \ q \-x+z-l I \ + G l {x-\,l, Z) 

1 \q\ x - z + 



r |g|*-*-3 + |g|-*+*+3 J 1 + 6?i(z + J, /, z) 



+ \q\~ x+z ){\q\^ + |g|~^-2r- 1 ) 



Ve 2 - 4(| 9 | 3i -*-3 + |g|-*+*+5)(|qf-*+3 + |g| -!B+ *-3) 



1 - 



a 



1 + - 1,1, z) 
[1-Gi(x,l,z)}. 



1 + Gi(.T + \,l,z) 



Hence we obtain bounds of the density one-point function 



(A.35) 
(A.36) 



-J^=) [1 - G~ 2 {x, I, z)\ < (n x ) c <(l- -A= ) [1 + Gt(x, I, z)\ (A.37) 



These Gf(x, I, z) consist of terms which decay exponentially for I — \x\, I — \z\ or \x — z\. 
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In the infinite- volume limit, these functions are asymptotically equal to 

lim Gf(x, I, z) = 0{p~ 2lx - zl ), (A.38) 

l— *oo 

for large \x— z\. Inequality (JA.35J) shows < 1, and hence the density one-point function 
does not vanish for all x G A. 



A. 4 Estimates of the Electron Two-Point Function 

We can represent the electron two-point functions in the following form 

(4 gCvT ) c =K'(x,y,l,z)—= 2g^ %= ~— , (A.39) 

' " y/l + \q\ 2 (*-*Wl + \<l\ 2iv - z) \ r 7 

where 



K'(x,y,l,z) = —> 



1 + 


'7 


2(1/-*) 




'7 


2(i-z) 



1 + 


( 7| 


2(w-z-±)} 


B{-l,x-l,z)B(y + l,l,z) 


1 _|_ g2M 1 


B{-l,l,z) 



(A.40) 



for x,y G A G . We can easily obtain K'(x,y,l,z) for other x,y E A, but they are com- 
plicated forms. Since the rescaled normalization function converges in the limit I — > oo, 
K'(x, y, /, z) exists in the infinite-volume limit. Apparently, K'(x, y, I, z) seems to diverge 
according to the factor ^/l + |g| 2 ^~ z ). But we can evaluate 

f |„|-2(i»-*)+l i 1 

igi-^ n |g| 



w=u+l 



1 + 


9 


2(w-*-£) 


1 + 1 


9 


2(uj-jz) 



and then the factor 



1 |q|-2(«J-z) _|_ 1 
u+1 11 

f |_|-2(i»-*)+2 ,1 1 i \ n \-2(v-z) 

<\qr {y ~ v) TT 11 , I \ =\q\- (y - v) + 9 , , (A.41) 

11 |g -2(u;-*) _|_ J l y l 1 _|_ |g|-2(j/-z)' V / 



1 + 


9 


2(lf-a) 




q 


2(x-z) 



1 + 


9 


|2(»-z-±) 


1 + 1 


9 


2(iu-z) 



9 


-2(y-v 


> + 


9 


2(v-z) 


1 + 


'7 


2(x- 





1 + 


'7 




1 + 


9 


-2(u-a) 


1+1 


9 


2(io-z) i 


1 + 


9 


-2(if-«) 



< \l 1 , I_I9.I%-^ I ll 1 . l-.l9r,„_^ I 1 . I_l-9f,,-,>1 (A. 42) 

is finite for any x and j/. Thus, we have found that K'(x, y, I, z) is finite for x, y G A Q even 
in the limit of I — > oo. We can show that K'(x,y,l, z) is finite also for other x, y G A. 
The expression (JA.39|) shows that the electron two-point function decays exponentially as 
y — x becomes large with a decay length 1/ logr because of the finiteness of K'(x, y, I, z). 
Thus, the electron two-point function (cj. a Cy,r) satisfies the requirement of the cluster 
property defined in subsection 5.5. 
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A. 5 Results of the Density Two-Point Function 



We can estimate the density two-point function as well as the one-point function. The 
density two-point function is written in 



A 4 B(-l, x - 1, z)B(x + 1, z)B(y + 1, I, z) 

T 2 B(-l,l,z) 



(A.43) 



for x, y £ A c . We can easily obtain representation for other x, y £ A, but it has a 
complicated form. The representations in terms of the rescaled normalization function 
give us good estimates. From the relation (|84j) . we obtain 

B(-l, I, z)=(l + \) B(-l, x-1, z)B(x + l,y-l, z)B{y + 1, /, z) 



l-f(x-\- z ) 



B{-1, x - 2, z)B(x + 1, y - 1, z)B(y + 1, /, z) 



+ 



l-f(x + \-z) 
l-f(y-l-z) 
l-f(y + \-z) 



B(-l, x-1, z)B(x + 2,y-l, z)B(y + 1, 1, z) 
B(-l, x-1, z)B(x + 1, y - 2, z)B(y + 1, /, z) 



B(-l,x- l,z)B(x + l,y - l,z)B(y + 2,l,z) 



+ 
+ 



1 - f(x - \ - z) 1 - f(y - \ - 











1 


- f(x - \ - z) 


1 


-f{y + \-z) 










1 


-/(z + i-z) 


1 


- f(y- \ - z) 










1 


-f(x + l-z) 


1 


-f(y+l-z) 











B(-l, x-2, z)B(x + l,y-2, z)B(y + 1, 1, z) 

B{-1, x-2, z)B(x + l,y-l, z)B(y + 2, 1, z) 

B(-l, x-1, z)B(x + 2,y-2, z)B(y + 1, 1, z) 

B(-l, x-1, z)B{x + 2,y-l, z)B(y + 2, 1, z). 

(A.44) 



From this relation, we evaluate the density two-point function 

' A 4 



\im(n x n y ) ( = < 

l— >oo 



2 _ G 3 (x,y,z) 
A 2 ^ 2 



if x, y £ A 
G4,(x,y,z) ilx,yeA' 

Ve — 4/ 
A 2 A 2 



e 2 -4 
1 



(A.45) 



1 - 



G 5 (x, y, z) otherwise 



Functions G3, G4 and G5 are definite functions which approach to 1 exponentially as 
\y — x\, \y — z\ and \x — z\ become large. 
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